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Abstract 

We prove the existence of a continuous family of positive and generally non-monotone 
travelling fronts for delayed reaction-diffusion equations ut{t, x) = Au{t, x)—u{t, x)+ 
g{u{t — h,x)) (*), when g £ C^(R+,M+) has exactly two fixed points: xi = and 
X2 = K > 0. Recently, non-monotonic waves were observed in numerical simula- 
tions by various authors. Here, for a wide range of parameters, we explain why 
such waves appear naturally as the delay h increases. For the case of g with nega- 
tive Schwarzian, our conditions are rather optimal; we observe that the well known 
Mackey-Glass type equations with diffusion fall within this subclass of (*). As an 
example, we consider the diffusive Nicholson's blowflies equation. 
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1 Introduction 



In this paper, we study the existence of positive non-monotone traveUing waves 
for a family of delayed reaction-diffusion equations which includes, as a par- 
ticular case, the diffusive Nicholson's blowflies equation 

Nt{t, x) = dAN{t, x) - 6N{t, x) + pN{t - h, x)e-^^^^-^'''\ (1) 

t e R, x e W^. This problem was suggested in [9,10,12,20,29], where numerical 
simulations indicated a loss of monotonicity of the wave profile caused by the 
delay. Eq. (1) was introduced in [30] and it generalizes the famous Nicholson's 
blowflies equation 

N'{t) = -5N{t) + pN{t - /i)e-''^^*-'^), (2) 

intensively studied for the last decade (e.g. see our hst of references). After 
a linear rescaling of both variables and we can assume that 5 — h — 1. 
Eq. (1) takes into account spatial distribution of the species, and the men- 
tioned problems reflect the interest in understanding the spatial spread of 
the growing population [13]. Relevant biological discussion can be found in 
[1,9,11,12,20,29], where various modifications of (1) were proposed and stud- 
ied. Here, however, we will concentrate mainly on the mathematical aspects 
of the dynamics in (1). For the sake of simplicity, we will consider the case 
of a single discrete delay, but extensions for more general functionals (which 
additionally can take into account non local space effects) are possible (cf. 
[6,21,22]). Since the biological interpretation of N is the size of an adult pop- 
ulation, we will consider only non-negative solutions for (1) and for other 
population models. Actually, our approach allows us to study a more general 
family of scalar reaction-diffusion equations 

ut{t, x) = dAu{t, x) - u{t, x) + g{u{t - h, x)), u(t, a;) > 0, a; e M™, (3) 
related to the Mackey-Glass type delay differential equations, 

u'{t) = -u{t) + g{u{t-h)), u>Q, (4) 

with exactly two non-negative equilibria Uiit) =0, U2{t) = K > Q (so that 
g{K) = K, g{0) = 0). In particular, with g{u) = pu/{l +-u") in (4), we obtain 
the equation proposed in 1977 by Mackey and Glass, to model hematopoiesis 
(blood cell production). The non-linearity g is called the birth function and 
thus it is non-negative, and generally non-monotone and bounded. Due to 
these properties of g and the simple form of dependence on the delay in (3), 
the Cauchy problem 

u{s, x) = C{s, x), s e [-h, 0], xe W^, (5) 
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for Eq. (4) has a unique eventually positive global solution for every C 7^ 
taken from an appropriately chosen functional space (e.g. see [27]). 

Recently, the existence of travelling fronts connecting the trivial and posi- 
tive steady states in (1) was studied in [31] and [6] (see also [9,23] for other 
methods which eventually can be applied to analyze this problem). In [31], 
the authors use a monotone iteration procedure coupled with the method 
of upper and lower solutions. This approach (proposed in [32]) works well if 
1 < < e, since in this case the function g is increasing on [0, 1/6] D [0, -ft'], 
thus —50(0) + pg{(l){—h)) satisfies the quasimonotonicity condition in 
[32]. This allows one to establish the existence of monotone wave front solu- 
tions N(t,x) = (f){ct + u • x,c) for every p/S E (l,e] and c > 2^yp — 5 (cf. 
[25,31]). Moreover, as it was proved in [25], every solution of (1), (5) with 
p/6 G (1, e] converges exponentially to some travelling wave provided that C is 
sufficiently close (in a weighted norm) to this wave at the very beginning 
of the propagation. For the case p/S > e, clearly g is not monotone on [0, K], 
and Wu and Zou's method [32] is no longer applicable. In [6], the Lyapunov- 
Schmidt reduction was used to study systems of delayed reaction-diffusion 
equations with non-local response. We observe that Eqns. (1) and (3) fit into 
the framework developed in [6]. This approach requires a detailed analysis of 
an associated Fredholm operator and the existence of heteroclinic solutions 
of (4) (in [6], the latter was established with the use of the monotone semi- 
fiows approach developed by H. Smith and H. Thieme [26]). As a result, it 
was proved in [6] that, even when p/S > e, (1) possesses a family of travelling 
waves if Sh e (0, r*) for some r* < 1 (which is given exphcitly). The rather 
restrictive condition Sh < r* < 1 from [6] was considerably weakened in [8] 
by invoking a Schauder's fixed-point argument to find heteroclinic solutions of 
(4). Unfortunately, the main results of [6,8] do not answer the question about 
the existence (and shape) of positive travelling fronts of (1) or (3). We recall 
here that only non-negative solutions to (3) are biologically meaningful. 

In this paper, inspired by [6,31,32], for a broad family of nonlinearities g (which 
includes Eq. (1) with S — 1), we prove that Eq. (3) has a continuous family of 
positive travelling wave fronts u{t, x) — (f){ct -\-v-x,c), indexed by the speed 
parameter c > 0, provided that 

e-'>-rin^^, Vilg'iK), (6) 

and c is sufficiently large: c > c^{h,g\Qi),g'{K)). Furthermore, we show that 
these fronts generally are not monotone: in fact, they can oscillate infinitely 
about the positive steady state. On the other hand, for large negative values of 
s, the wave profile 0(s, c) is asymptotically equivalent to an increasing expo- 
nential function. Condition (6) assures the global attractivity of the positive 
equilibrium of (4), which is required by our approach. It should be noted that 
this condition is rather satisfactory in the sense that (6) determines a domain 
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of parameters approximating very well the maximal region of local stability 
for the positive steady state in (4) or (2) (cf. [22]). 

Before announcing the main results of the present work, we state our basic 
hypothesis: 

(H) Eq. (4) has exactly two steady states ui(t) = and U2{t) = K > 0, the 
second equilibrium being exponentially asymptotically stable and the first 
one being hyperbohc. Furthermore, g e C-'^(]R_,_, ]R_|_) and is C^- smooth in 
some vicinity of the equilibria, with p := g'{0) > 1. The latter implies that 
the solution ui — of (4) is unstable for all /i > 0. 

In the sequel, Ai(c) denotes the minimal positive root of the characteristic 
equation {z/cy — z — 1 + pexp{—zh) — for sufficiently large c, and A the 
unique positive root of the equation —z — 1 + pexp{—zh) = 0, where p > 1. 
As shown later, limc^oo Ai(c) — X. Now we are ready to state our main result: 

Theorem 1 Assume (H). If the positive equilibrium K of Eq. (4) is glob- 
ally attracting, then there is c^, > such that, for each v G M™, \v\ = 1, 
equation (3) has a continuous family of positive travelling waves u{t,x) = 
(f){ct + u ■ x,c), c > c*. Furthermore, for some sq = so(c) G M, we have 
0(s — So, c) = exp(Ai(c)s) + 0(exp(2As)) as s — > — oo, so that 4)'{s — sq, c) = 
Ai(c) exp(Ai(c)s) + 0(exp(2As)) > on some semi-axis {—oo,z]. Finally, if 
g'{K)he^^^ < — 1 then the travelling profile (j){t) oscillates about K on every 
interval [2;, +oo). 

In order to apply Theorem 1, one needs to find sufficient conditions to ensure 
the global attractivity of the positive equilibrium of (4). Some results in this 
direction were found in [21,22] for a family of nonlinearities having negative 
Schwarz derivative (or, more generally, satisfying a generalized Yorke condition 
[7,21,22]). In particular, [21, Corollary 2.3] imphes the following useful version 
of Theorem 1: 

Corollary 2 Assume (H) and (6). In addition, suppose that g G C"^(M+,M_|_) 
has only one critical point xm (maximum) and that the Schwarz derivative 
{Sg){x) = g"'{x){g'{x))-^ - (3/2) {g" {x){g' {x))-^f is negative for all x > Q, 
X 7^ Xm- Then all conclusions of Theorem 1 hold true. 

Notice that Corollary 2 applies to both the Nicholson's blowflies equation and 
the Mackey-Glass equation with non-monotone nonhnearity, see [21]. 

To prove our main results, we need a detailed analysis of heteroclinic solutions 
of (4). This study is presented in Section 2, and is crucial for the selection of an 
appropriate functional space where a Lyapunov-Schmidt reduction is realized. 
The existence of positive travelling waves is proven in the third section. The 
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main result of Section 3 is given in Theorem 14, which is essentially Theorem 

1 without its non-monotonicity statement. Finally, in the last short section, 
we show that these waves have non-monotonic profiles when the delay is over 
some critical value. 

2 Heteroclinic solutions of scaleir delay differential equations 

In this section, we study the existence and properties of heteroclinic solutions 
to the scalar functional equation 

x'{t) = -x{t)+f{xt), x>0, (7) 

where / : C([— /i, 0], M+) R+ is a continuous functional which takes closed 
bounded sets into bounded subsets of R+. Here C([— 0], IR+) is the metric 
space equipped with the norm |0| = maXsg[_/i^o] Throughout this pa- 

per, we suppose that Eq. (7) has exactly two steady states Xi{t) = and 
X2{t) = K, the second equilibrium being asymptotically stable and glob- 
ally attractive. Thus, if (7) has a heteroclinic solution it must satisfy 

lj;(-Oo) = 0, = K. 

We start by proving a general existence result which is valid for the abstract 
setting of dynamical systems. Let : X ^ X he a. continuous semidynam- 
ical system defined in a complete metric space {X,d). First, we mention the 
following fact (see e.g. [16, p. 36]): 

Lemma 3 Suppose that 99 : M — > X, 99(0) = x is a complete orbit of S^. 
If the closure {(p{s), s < p} is compact for some p e R, then the a— limit 
set a{ip) = ng<o{(f{s), s < q} of 1^ is nonempty, compact and invariant (this 
means that for every z G q;(<^) there exists at least one full trajectory ip with 
i;{R)C a{ip), i^{Q) = z). 

For every A C X and h > 0, let A{h) C A denote the set of right endpoints of 
all orbit segments S^^'^^z — {S'^z : u E [0,h]} which are completely contained 
in A: 

A{h) = {x e A:x^ S'^zand S^'^^'^h C A, for some z e A}. 

Next statement shows clearly how to relate the global attractivity property 
of the positive equilibrium of (4) to the problem concerning the existence of 
travelling fronts for (3): 

Lemma 4 Assume that A{h) is either empty or pre- compact, for all bounded 
sets A and some h > 0. Suppose that there exist two disjoint compact invariant 
subsets Ki,K2 of X such that d{S*x,K2) — >■ as t — >• -|-oo for every x G 
X \Ki. If the set — {x : d{x, Ki) = e} is not empty for every sufficiently 
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small £ > 0, then there exists at least one complete orbit ip with a{ip) C Ki 
and u!{ip) C K2. 



PROOF. Let p = d{Ki,K2) and, for every n > 2/p, take some G 
Fi/n- Due to the compactness of Ki, we can assume that Xn ^ z for some 
z e Ki. In consequence, if i„ > is the minimal real number such that 
d{S*"Xn,Ki) = p/2, then limt„ = +00. Set Wn — S^^-Xn- Due to the com- 
pactness condition imposed on S**, we can suppose that limw„ = w. Let 
now ipniu) = S^~^^"Xn, u > —tn- Wc havc S'^ipnif) = i^ni.O' + t) for every 
a > 0,t > —tn- Since, for every integer m > the sequence ipn^—m) has 
a convergent subsequence (say, ipuj {—'itT') b), we can assume that ipnit) 
converges uniformly on [— m, 0] to ilj{t) = S*~^"^b. Moreover, we have that 
il^{0) = w and S"'il){t) = il){a + t) for all a > 0,t > — m. In this way, tak- 
ing m = 1,2,3,..., we can use ipn{u) to construct a continuous function 
: M — > X, such that S"'ip{t) = ip{a + t) for every a > 0, t G M. Such ip defines 
the complete orbit we are looking for. Since '0(M_) is a subset of the bounded 
set B — {z : d{z,Ki) < p/2}, we conclude that '?/'(R_) is pre-compact. Fur- 
thermore, because of d{ip{W_), K2) > p/2, we have d{a{ilj),K2) > p/2 > 0. 
This means that a{ip) C Ki so that limt^_^ d{ip{t), Ki) = 0. □ 

A direct application of Lemma 4 to Eq. (7) gives the following 

Theorem 5 Let f : C([— /i, 0], M+) R+ be a continuous functional which 
takes closed bounded sets into bounded subsets of M+. Assume further that 

every non-negative solution of (7) admits a unique extension on the right semi- 
axis. // /(O) = 0, f{K) = K {K > 0) and X2(t) = K attracts every solution 
of (7) with non-negative and nontrivial initial function, then there exists a 
positive complete solution iplt) to (7) such that ip{—oo)=0, ip{-\-oo) — K . 

With some additional conditions on /, we can say more about such an orbit 

Lemma 6 Assume f{(f)) — g{(f){—h)) for some g G C(R+,R+) (so that g{0) — 

0, g{K) = K). Assume fkatXimmix-^Qj^ g[x)/x > 1, and letpi,p2 be such that 
1 < Pi < \iminfx->o+ gi^) /x < limsup^^g_|_5'(a;)/a; < p2. Let Aj be the unique 
positive real root of the equation z = —1 -\- Piexp{—zh) (i = 1,2), so we have 
< Ai < A2. Then for every heteroclinic solution ifjit) of the equation 

x'{t) = -x{t) + g{x{t - h)) (8) 

there exist r = T{ip) < 0, Cj = Ci{ip) > such that 

Ciexp(A2i) < i/^it) < C2exp(Aii), t < r. 
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PROOF. Choose S > sufficiently small such that pix < g{x) < p2X for all 
X G [0, 6), and let r be such that ip{t) < S for all t < t. We claim that, for 

every s <t 

/ / \ / / \ exp(— /i) exp(— /i) 
'0(^m) = mm ■i/y(s + M)> max %p[s + u) — Wv-m)- (9) 

u6[-/i,0] P2 M6[-/i,0] P2 



Indeed, if tm > tM, then, by the variation of constants formula, 

■ip{tm) = ^{tu) exp{tM-tm)+ / exp{-{t^-u))g{ilj{u-h))du > ^(^m) exp{-h). 



tM 

Finally, suppose that tM — h < tm < so that ip'{tM) > 0. Then (9) holds 
since 



ipitn) < g{ip{tM - h)) < p2ip{tM - h); ip{tm) = ip{tM - h) exp(tM -tm-h)^ 
j exp(-(t^ - u))g{'ilj{u - h))du > il^itM - h) exp(-/i). 

tM—h 

Next, for every s < t and u e [— /i, 0], we have that 

ii{tm) exp(AiM) < V'(s+m) < '0(tM) exp(A2(M+/i)) < P2i>{tm) exp{X2{u+h)+h) . 

Prom the inequalities above and since additionally piX < g{x) < p2X for all 
X e [0, S), then ior s + h < r and e [— /i, 0] we have 

s+h+u 

tlj{s + h + u)^ ^{s)e-^^+''^ + e-(^+'^+") j e^git/jia - h))da 

s 

s+h+u 

< z/;(s)e-(''+") + e-(^+'^+"V2^(tM) / e'^+^^Ca-.) 

and 
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By repeating the above procedure over intervals of length h, the step by step 
method implies that, for all —h < u < t — s, 

'^{tm) exp(Aiu) < ^'(s + u) < P2'^{tm) exp(A2(w + /l) + /l). 
In particular, 
ipis) 

exp(Ai(r-s)-/i) < V'(^m) exp(Ai(r-s)) < ^/'(r) < p2^(tm) exp(A2(r-s+/i)+/i). 

Thus, for every s < r, 

P2^iP{t) exp{X2{-T-h+s)-h) < ip{tm) < V'(s) < P2ij{T) exp(Ai(-r+s)+/i). □ 



In what follows, we shall assume that g G C"'^(]R_|_, M+), g'{0+) — p > 1, and 
use several times the following simple assertion. 

Lemma 7 Suppose that p > 1 and h > 0. Then the characteristic equation 
z = —1 + pexp{—zh) (10) 

has only one real root < \ < p — 1. Moreover, all roots A, A^, j = 2, 3, . . . of 
(10) are simple and we can enumerate them in such a way that 

. . . < 5RA3 = 5RA2 < A. 



PROOF. The last inequality follows from "RXj < -1 +pexp(-/i3?Aj), j > I. 
□ 

Lemma 8 Suppose that 5''(0+) = p > 1 and that ip is a heteroclinic solution 
to (8). Let A he the positive root of (10). If there exists g"{^+) G M, then, for 
each 5 > and some to e we have that ip{t — to) — exp{Xt) + 0(exp((2A — 
S)t)) at t -00, so that ip'{t - to) = Aexp(At) + 0(exp((2A - S)t)) > 
on some semi-axis (— oo,T]. Moreover, if there exists g" : [0, e) R and is 
bounded for some e > 0, then ip{t) is unique up to a shift in t. 

PROOF. Since ^"(0+) is finite, g{x) = g'iO+)x + ©(x^) ^ px + 0{x^) as 
X — 0. From Lemma 6, given S > small, for Pi = p — S/2 we have ilj{t) — 
0{cxp{\st)) at —00, where A^ is the unique positive root oi z = —l + {p — 
6/2) exp{—zh). It is easy to see that A5 > A — 6/2. In fact, let W{z) :~ 
z + l — {p — 5/2) exp{—zh) . For 2; e R, we have W{z) < if and only if 2; < A5. 
On the other hand, W{X-5/2) = pexp(-A/i)[l -exp(5/i/2)] +5/2[exp(-(A- 
6/2)h) - 1] < 0. Hence, 

^'(t) = -h) + 0{^l;'{t - h)), (11) 
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where 0{ip'^{t — h)) = 0(exp((2A — 5)t)) as i — > — oo. Now, consider the hnear 
inhomogeneous delay differential equation 

x'{t) = -x{t) +px{t -h) + 0(exp((2A - S)t)). (12) 

The change of variables x{t) = y{t) exp((2A — 5)t) transforms it into 

y'{t) = -(1 + 2A - 5)y{t) +pexp(-(2A - 5)h)y{t - h) + 0{1). (13) 

The spectra o-{y), a{x) of the hnear parts of Eqns. (13) and (12) are related by 
a{y) = cr[x) —2X + S, therefore the linear part of (13) has not pure imaginary 
eigenvalues for all sufficiently small 6 > (equivalently, the linearization of 
Eq. (13) about zero is hyperbolic). In this case, (13) has a bounded solution 
yb{t) = 0(1) at t = —oo (e.g. see [4, Lemma 3.2, p. 246] or [17, Section 
10.1]. Note that Eqns. (12) and (13) are not autonomous. Nevertheless, the 
results for autonomous equations near hyperbolic equilibria in [4,17] are valid 
in this setting, since the linearized equation near zero for Eq. (13) has an 
exponential dichotomy, cf. [17, p. 312]). Thus Eq. (12) has a solution xi){t) = 
yb{t) exp((2A — 6)t) = 0(exp((2A — S)t)). In consequence, z{t) = ■ip{t) — Xh{t) 
solves the linear homogenous equation x'{t) — —x{t)+px{t — h) and is bounded 
at i — > — oo. This is possible if and only if 

N 

z{t) = Cexp(Ai) + J2 ^3 exp(Aji), 

where A > 0, Aj e C, j = 1, . . . , is a finite set of roots having non- negative 
real parts of the characteristic equation (10). Notice that C e R, Cj e C and 
A > 3fJAj (see Lemma 7). In this way 

N 

ip{t) = Cexp(Ai) + Cj exp(Aji) + 0(exp((2A - 5)t)). 
i=i 

On the other hand, from Lemma 6 we know that ip{t) = 0(exp((A — S/2)t)). 
Since A > 3?Aj, this implies immediately that all Cj = 0, C > and that 
V^(i) = Cexp(At) + 0(exp((2A - 5)t)). By (11), 

jjj'(t) = CAexp(Ai) + 0(exp((2A - 5)t)) > 0. 

Observe also that ^{t) = '^(t-A-MnC) = exp(At) + 0(exp((2A-(5)t)) defines 
another heteroclinic solution of (7). 

Finally, suppose that ii{t),u{t) are two heteroclinic solutions to (7) such that 
H{t) = exp(At) + 0(exp((2A - 5)t)), u{t) = exp(At) + 0(exp((2A - 5)t)). 
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Applying the Lagrange mean value theorem twice, we get g{x) — g{y) = p{x — 
y) + {x — y)0{x + y) for y close to 0. Since (T{t) — ii{t) — — 0(exp((2A — 
S)t)) we obtain that 

g{fx{t-h))-g{iy{t-h)) = a{t-h){p+0{exp{Xt))) = pa{t-h)+0{exp{{3X-d)t)). 
Therefore a{t) satisfies 

x'{t) = -x{t) + px{t - h) + 0{exp{{3X - S)t)), (14) 

from which, applying the same procedure as above, we deduce that a{t) — 

/j{t) — v{t) = 0(exp((3A — S)t)). In this way, we can show that a(t) ~ 
0(exp((A;A — S)t)) for every integer k > 2. This leads us to the conclusion 
that a has superexponential decay at t = — oo (equivalently, a is a small solu- 
tion at t = — oo, see [4]). We will finalize our proof showing that only the trivial 
solution of the linear asymptotically autonomous homogeneous equation 

x'{t)^-x{t)+p{t)x{t-h), p(-oo)=p>l, (15) 

can have superexponential decay at t = — cxo (notice that a{t) satisfies (15) 
with p(t) = p + 0(cxp(At))). Indeed, if x{t) > on some semi-axis (— oo,^], 
then we can repeat the arguments in the proof of Lemma 6 to find an ex- 
ponential lower bound for x{t), in contradiction to our assumption of super- 
exponential decay of x{t). Consider now the case of x{t) oscillatory on every 
semi-axis (— oo,^], and take zq such that p{t) < p -|- 1 for all t G (— oo,2;o]. 
Let ti G (— oo,zo) be a point of the global maximum of |a;(t)|: we can as- 
sume that x{ti) = M > 0, x'{ti) > 0. Then x{ti - h) > M/{p + 1), 
so that |a;(t2)| = ^^^t<ti-h\x{t)\ > (p + 1)^"*^ niaxi<ija;(t)|. Analogously, 
maxt<ti_2h\xit)\ > \x{t3)\ = maxt<t2-h\xit)\ > {p + maxt<t2 1^(^)1 = 
(p-M)~^maxt<tj_,j \x{t)\ > (p-M)~^maxt<tj \x{t)\. Thus 

max \x(t)\ > (p+ 1)"'^ max |a;(t)| 
t<ti-kh' ^ ^1 — ^-^ ' t<ti ' ^ " 

so that x{t) can not decay superexponentially as i — > — oo. □ 

Now, assume (H) and the global attractivity of a;2 = -f^ for Eq. (8), and then 
take A > satisfying (10) and the unique (up to a shift in time) heteroclinic 
solution i/j described in Lemma 8. Let A* G (0, A) be sufficiently close to A 
and such that the equation y'{t) = —(1 + A*)y(t) + pexp(— A*/i)y(t — h) is 
hyperbolic. Note that this latter equation is obtained by effecting the change 
of variables a;(t) = cxp{X^t)y{t) to the linear equation a;' (t) = —x{t)+px{t—h). 
For a fixed // > 0, we will consider the seminorms = supjj^ 1^(^)1) ll^ll^ — 
sup]K_ e~^^\x{s)\, ll^llju = max{||a;||"^, and the following Banach spaces: 

Cju(M) = {x G C(M, M) : < oo and x(-|-oo) exists and is finite}. 
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C^,xM) = e Ca,(M) : J x{s)i:'{s)ds = 0}, 

— oo 

equipped with the norms and ||a;||A, respectively (in order to simpUfy the 
notation, we shall often write instead of etc). Notice that, due to 

Lemma 8, we have i/j' G Ca. {^) \ (M). We shall also need the following 
integral operator 

t 

U : C^,A,(M) ^ Ca,(M); {Mx){t) = j e-^'-'h{s)x{s - h)ds, 

— oo 

where q{s) = g'{ip{s — h)) with g(— oo) = g'{^+) = p > 1, g(+oo) = g'{K). 
Observe that M is well defined, since {Nx){+oo) = g'{K)x{+oo) and, for 
t < h, 

mx){t)\ = / e-^'-^^\q{s)\\\x\\^^e^*^^-^^ds < tK^El^lMAe^*^'-'^\ 



Lemma 9 //(H) is assumed, then I —J\f : C^,a*(IR) C\^{M.) is an isomor- 
phism of Banach spaces. 

PROOF. We first prove that Ker(/-A/') = 0. Indeed, if y e Ker(/-A/') and 
y 7^ 0, then 

t 

J e-'^'-'k{s)y{s-h)ds = y{t). 

— oo 

Therefore y is a bounded solution of the linear delay differential equation 
y'{t) = -y{t) + q{t)y{t-h). (16) 

Since g'{x) = p + 0{x) at x = and ^'(O — 0{exp{Xt)) at t — — oo, we 
conclude that 

q(t) = g'{'ilj{t -h))^p + (9(exp(Ai)), t -oo. 

Thus y{t) can be viewed as a bounded solution of the inhomogeneous equation 

x'{t) = -x{t) +px{t -h) + 0(exp(2At)), t -oo. 

Since y{t) = 0(exp(A*t)) at — oo, with A* < A close to A, the procedure which 
has been used before to prove the uniqueness of the heteroclinic ipit) allows us 
to conclude that |/(t) = Cexp(At)+0(exp(2A*t)) andthat dimKer(/— A/") = 1. 
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On the other hand, wc know that ijj'it) ^ satisfies (16). Thus we must have 
y(t) = cip'it) ^ C^^A. (IR), c 7^ constant, a contradiction. Therefore y{t) = 
and Ker(/-7\A) = 0. 

We now estabhsh that / — A/" is an epimorphism. Take some d e (R) and 
consider the following integral equation 

t 

x{t)- j e-^*-'\{s)x{s-h)ds^d{t). 

—oo 

If we set z{t) — x{t) — d{t), this equation is transformed into 

t 

z(t) - J e-^^-'^q(s)(z(s - h) + d(s - h))ds = 0. 

— oo 

Hence we have to prove the existence of at least one Cx^ (R)-solution of the 
equation 

z'{t) = -z(t) + q{t)z{t -h) + q{t)d{t - h). (17) 

First, notice that all solutions of (17) are bounded on the positive semi- axis 
R+ due to the boundedness of q{t)d{t — h) and the exponential stability of the 
homogeneous a;-hmit equation z'{t) — —z{t) + g'{K)z{t — h). Here we use the 
persistence of exponential stability under small bounded perturbations (e.g. 
see [3, Section 5.2] or [5, Chapter VI (9c)]) and the fact that g(+oo) = g'{K). 
Furthermore, since every solution z of (17) satisfies z'{t) = —z{t) + g'{K)z{t — 
h) + g\K)d{+oo) + e{t) with e(+oo) = 0, we get ^(+oo) = d{+oo)g\K){l - 
g'{K))~^. Next, by effecting the change of variables z{t) — exp(A*i)y(i) to Eq. 
(17), we get a linear inhomogeneous equation of the form 

y'{t) = -(1 + K)y{t) + [pexp(-A,/i) + ei{t)]y{t - h) + e2(t), (18) 

where ei{—oo) = and e2(t) = 0(1) at t = —oo. Since the a-limit equation 
y'{t) — —(1 + X*)y{t) + pexp(—X^,h)y{t — h) to the homogeneous part of (18) 
is hyperbolic, due to the above mentioned persistence of the property of ex- 
ponential dichotomy, we again conclude that Eq. (18) also has an exponential 
dichotomy on ]R_. Thus (18) has a solution y* which is bounded on R_ so that 
z*{t) = exp{Xj)y*{t) = 0{exp{Xj)), t -oo, is a Ca, (M)-solution of Eq. 
(17). Now, it is evident that w{t) = z*{t) - Cijj'it) = 0(exp(A,t)) solves (17) 
for each C e M. In consequence, x{t) = d{t)+z*{t)-Cdip'{t) = {{I-J\f)-^d){t), 
if we take C, = + z*{sW{s))ds{f_^{i;'{s)rds)-\ □ 

Remark 10 For 6 > small, consider I — Ni : C2a-5(R) — >■ C2\-&{^), where 
Ml is defined by {Mix){t) = p e''^*~^^x{s — h)ds (recall here the discussion 
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after formula (13)). Replacing J\f hy Mi in the proof of Lemma 9, we establish 
similarly that I — Mi is an isomorphism of the Banach space C2x-s{^) onto 
itself. Since the linear equation x'{t) = —x{t) + px{t — h) is hyperbolic, this 
situation is actually simpler than the one considered in Lemma 9. 

3 Existence of a continuous family of positive travelling waves 

In this section, we are looking for travelling waves for (3), that is, solutions 

u{x,t) = ({){eh' ■ X + 1), z/ G M"*, = 1, where c= 1/e is the wave speed, 
connecting the two equihbria of (3). We will suppose that e is sufficiently small. 
This leads us to the question about the existence of heteroclinic solutions to 
the singularly perturbed equation 

e'^x"{t) - x'{t) - x{t) + g{x{t -h)) ^0, te M, (19) 



with x(— oo) = 0, x{+oo) = K. Being a bounded function, each travelling 
wave should satisfy the following integral equation 



x(i) = — — I y e^+^g{x{s-h))ds+ J e ^ g{x{s - h))ds \ , (20) 



where a{£) = \/T+~Ae^. For solutions in Ca^(M) with A* G (0, A) close to A, 
this equation can be written in the shorter form 

X - (J, o g)x = 0, (21) 
where I„g : Ca.(M) ^ Ca,(M) are defined by 

/ t +00 \ 

I I r -2(t-s) f (l+CT(e))(t-s) \ 

{X^x){t) = ——\ / e x(s - /i)(is + / e ^ x{s-h)ds\, 

"^'^ vL { J 

and {Qx){t) = g{x{t)) is the Nemitski operator. (For the sake of simplicity, we 
write Xe, Q instead of le^x^, Q\^). Throughout all this section, we will suppose 
that the C^-smooth function g is defined and bounded on the whole real axis 
R. Clearly, this assumption does not restrict the generality of our framework, 
since it suffices to take any smooth and bounded extension on M_ of the 
nonlinearity g described in (H). Notice that, since there exists finite ^''(0), 
we have g{x) = xj{x) for a bounded 7 G C(R). Set 70 = supjgjg |7(a;)|. As 
it can be easily checked, \\Gx\\ < 7o||a^|| so that actually Q is well-defined. 
Furthermore, we have the following lemma: 

Lemma 11 Assume that g G C^(M). Then Q is Frechet continuously diff'er- 
entiable on Ca. (M) with differential Q\xq) : y{-) — > g\xQ{-))y{-). 
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PROOF. By the Taylor formula, g{v) — g{vQ) — g'{vo){v — vq) — o{v — vq), 
I'd e M. Fix some Xq G Ca»(IR), then we have 

WQx - gxo - g\xo{-)){x - xo)\\ = o{\\x - xoW), X e Ca.(R). 

Clearly, it holds that = ||5f'(x(-))M(-) || < sup^gj^ [(^'(^(t))! Since 

functions in Ca. (M) are bounded and g' is uniformly continuous on bounded 
sets of R, for any given S > there is cr > such that for ||a; — a;o|| < cr we 
have \\g'{x)-g'{xo)\\<S. □ 

Now, we consider the integral operators : Cy^(M) C'^(R) defined as 

+00 t 
= / e ^ x{s-h)ds, {I~x){t) ^ / e~^+^ x{s - h)ds. 
t —00 



Lemma 12 Set X — Xq and Xq =0. If e ^ 0+, then X^ ^ X in the operator 
norm. Moreover, both operator families X^ : [0,1/ y/JI) >C(C^(]R), C^(]R)) 
are continuous in the operator norm. 

PROOF. We prove only that {{X^ — ^\\ ^ as £ ^ 0, the proof of the 
continuous dependence of X^ on e being completely analogous. 

We first establish that X+ — > uniformly as £ — > 0. In fact, for all t e M, we 
obtain 

+00 

\{X+x){t)\< J e'-^\x{s-h)\ds<e^x\\. 
t 

Furthermore, since \x{t)\ < \\x\\ exp(/xt) for all i e R, for £^ < we have 

t-s \\x\\ 
',~^\x(s — hj\ds < — 7^ — 

Hence, for £^ < (1 — 0.5e~^'^)//x, we obtain that < 2£^||a;||. 

Next, we prove that — > X uniformly as £ — > 0. We have 



-TOO 11 11 

{X+x){t)\ < j e^\x{s-h)\ds< X_ g""^'"^^- 



—00 

Thus, for t < h, we obtain that \{{X~ — X)x){t)\ < 



(2 + (a(£) + l)//)(l + /.)' 
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and, for all t, 



\{{i--i)xm\< I e-(*-)( 



'^(e)-i 



its) 



— l)\\x\\ds — \\x 




2 



— oo 



Thus \\X^ — X\\ < 0.5{a{e) — 1), and the proof of the lemma is complete. □ 

To prove the main result of this section, stated below as Theorem 14, we will 
make use of the following proposition: 

Lemma 13 Let {za, o. e A}, where NU {oo} C A, denote the (countable) set 
of roots to the equation 



If p > 1, > 0, e G (0, l/(2-y/p — 1)) then (22) has exactly two real roots 
1 ^oo{^) such that 



Moreover: (i) there exists an interval O = (0, a(p, h)) such that, for every e e 
O, all roots Xa{£),o. e A of (22) are simple and the functions Aq : C — ^ C are 
continuous; (ii) we can enumerate Xj{e),j e N, in such a way that there exists 
lim^^Q^ Xj{e) = \j for each j G N, where \j G C are the roots of (10), with 
Ai = A; (Hi) for all sufficiently small e, every vertical strip ^ < "Rz < 2{p — 1) 
contains only a finite set of m{^) roots (if ^ ^ {KA^, j G N}, then m(^) does 
not depend on e) Xi{s), . . . , Xm(^){6) to (22), while the half-plane > 2(p— 1) 
contains only the root Aoo(£). 

PROOF. The existence of real roots Xi{e),Xoo{B) satisfying A < Xi{e) < 
Xoo{e) is obvious when e G (0, 0.5/ ^/p — 1). On the other hand, if 2:0 > is 
a real root of (22), then e'^Zq — 2^0 — 1 < 0, e^Zq — 2^0 — 1 + p > 0. Hence 
^0 < (1 + Vl + 4£^) / (2£^) < + 1, from which it can be checked easily that 

1 + ^1 - 4(p - 1)^2 i_2(p-l)£2 1 - ^1 - 4(p - 1)^2 

Aoo > > , Ai < — < 2(^ 

We also notice that every multiple root Zq has to satisfy the system 

£^2^0 — zo — 1 + pexp{—zoh) — 0, 2£^2;o — 1 — phexp{—zoh) — 0, (23) 

which implies 



e'^z'^ - z -1 +pexp{-zh) = 0. 



(22) 



< A < Ai(£) < 2(p- 1) < E 



-2 



2{p - 1) < Aoo(£) < + 1. 



{e'^Zq — Zo — l)h + 2e^zq — 1 = 0, pexp(— 2^0/1) 



2 + 2:0 
2 + hzQ' 



(24) 
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The first equation of (24) implies that zo is real while the second equation of 
(23) says that Zq > 0. Since Zq is positive, from the first equation of (24) we 
obtain 0.5s~^ < Zq (we recall that e^Zq—Zo—I < 0). Let Co{p, h) be the maximal 
positive root of the second equation of (24). If £ > is so small that O.Se"^ > 
Co{p,h), system (23) can not have any positive solution. In consequence, the 
second assertion of this lemma holds if we set a{p, h) — l/y^2Co(p, h). 

Finally, we prove that the half-plane > 2(p — 1) contains only the root 
Aoo(£^) of (22). For this, let us evaluate le^z"^ — z — 1\ on the boundary of 
some rectangle [2(p — l),h\ x [— c, c] C C, with 6, c being sufficiently large. For 
/^(s), v{e) the (real) roots of e^z"^ — 2; — 1 = 0, we have that 

\e^z^-z-l\ = s^\z-n{e)\\z-u{e)\ > e^\^z-ii{e)\\^z-u{e)\ = \e\^zf-^z-l 

Thus, for ^z — 2{p — 1), we obtain 

leV - ^ - 1| > 3fJ^ + 1 - s^^zf > p. 

If Uz > 2(£-2 + l), then 

Is'^z'^ - ^ - 1| > e\^zf -Uz-l>8p-3>p. 

Similarly, for \'tsz\ > p/s fixed, we get 

leV - z - 1| = e^\z - n{e)\\z - iy{e)\ > e^{^zf > p. 

Thus, by Rouchc's theorem, e'^z'^ — z — 1 +pexp{—zh) = and e'^z'^ — z — 1 = 
have the same number of roots in the half-plane ^z > 2(p — 1), that is exactly 
one root. 

Therefore, for aU with ^Xj e [^, 2(p - 1)] and e e (0, 0.25/ VP^), we get 

pe'^'' > |S>(£^A| - - 1)1 = |S>Aj||l - 2e^^Xj\ > \QXj\/2, 

so that I'^Xjl < 2pe^^^. Hence, applying Rouchc's theorem to the functions 
e'^z'^ — z — 1 + pexp{—zh) and —z — 1 + pexp{—zh) along an appropriate 
rectangle inside — 1, 2{p — 1)] x [—3pe~^^, 3pe~^'*] C C, we prove the last 
assertion of Lemma 13. □ 

Theorem 14 Assume (H) and that the positive equilibrium of Eq. (4) is 
globally attractive. Let ip be some heteroclinic orbit of Eq. (J,.): ip^—oo) = 0, 
il!{+oo) ~ K. Then, for every 5 > there is an interval £ — (— £0,^0) o,nd a 
continuous family of positive heteroclinic orbits tp^ : £ ^ Ca-5(M) of Eq. (19) 
satisfying the additional conditions = and J^^'ijj^{s)ip'{s)ds = 0.5'?/'^(0). 
Furthermore, for every e & £ \ {0} we have that ipi;(t — tg) = exp(Ai(£:)t) + 
0(exp(2At)) at t —00 for some to = to{e) G and that ■ip'^{t — to) ~ 
Ai(£) exp(Ai(£)t) + 0(exp(2At)) > on some semi-axis (—00,2;]. 
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PROOF. We represent the mentioned orbit tp of (4:) as tp — aip' + 4>o, where 



—oo 

For S > small, consider A* = A — 5. In virtue of Lemmas 9, 11 and 12, we 
can apply the implicit function theorem (e.g., see [2, pp. 36-37] or [28, p. 170]) 
to Eq. (21) written as F{(f),e) = 0, where F : C^,xJ{R) x R ^ Ca,(M), 

F{(f), s) = a^p' + 4)- (J, o g){a^p' + 0), and Xq = ^■ 

Observe that F(0o, 0) = and F^{(f)Q, 0) = I —M. In this way, we estabhsh the 
existence of an interval £ = (— eo,£^o)) ^ (0, l/{2^p — 1)) and a continuous 
family (f)^ : S ^ C^^x^{M.) of solutions to F{(j),e) = 0. Notice that ipQ = ijj, 
ips = aip' + (pe £ C'a, (-R) satisfy Eq. (21), so that, as it can be checked directly, 
'ijj£{+oo) = g{^jjs{+oo)). Thus tjjs{+oo) — K and -0^ satisfies all conclusions of 
the third sentence of the theorem, except its positivity, which is proved below. 

Assume now that is sufficiently small so that \i{e) < 0.5Xoo{s) for all 
e e £ \ {0}. Since g{x) = px + 0{x'^) as x 0, and since there exists a 
constant Ci > such that {ipeit)] < Ci exp(A*t), t <0, e e E, we get 

e'^':{t) - ij'M - ^P,{t) +pMt -h)^ ^S), (25) 



where ^',(-) = pVs(- - h) - g{M- - h)) e C2xM- Moreover, ||^r,||2A. < C2 
for some C2 > which does not depend on e. Now, C2A* (R)-solutions x^ to 

s^x"{t) - x'{t) - x{t) + px{t -h) = ^s(t), (26) 



are solutions to the equation (/ — pls)xs = —le^e- Due to Remark 10 and 
Lemma 12, for A* = A — 5 close to A the operator / — pX^ is invertible in 
C2A. (IR) for all sufficiently small e. Moreover, Lemma 12 implies that there 
exists a subinterval £-i<Z £ such that ||(/— pJ^)"^!! < C3 for all e e £1. Hence, 
we obtain \\xe\\ < \\{I -ple)-^le\\\\^e\\ < C4 for all e E £1. Therefore Eq. (26) 
has a bounded solution x^ such that |a;£(t)| < C4exp(2A*t), t < 0, e E £1. 
Consequently, Zs{t) = ip^{t) — Xs{t) solves the linear homogenous equation 

e^z'^t) - z'{t) - z{t) + pz{t - /i) = 0, teR 

and is bounded as t ^ — 00. This is possible if and only if 

N 

Ze{t) = ^eexp(Ai(£)i) + exp(Aoo (e)^) + Y^Cj^eQ'^v{\{.£)t), 

J=2 
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where \j{e) & C, j = 1, . . . , N and Xooi^) are the roots with non-negative real 
parts of the characteristic equation (22), A^, e M, Cj^^ € C, £ e £i \ {0}. In 
consequence, 

N 

xpe{t) = A^exp{Xi{e)t) + B^exp{Xoo{e)t) + Y,Cj,eexp{Xj{e)t) + Xe{t). 

It follows from Lemma 13 that '^Xj{e) < A* < A < Ai(£:) < 0.5Aoo(£), 
provided e is small (say, e e £2 ^i) and A* is sufficiently close to A. Since 
ipeit) — 0(exp(A*i)), this implies immediately that Cj^^ — and 

i/j.it) = A, exp(Ai(£)t) + B, exp(Aoo(£)t) + x,{t), t G M, £ G £2 \ {0}. 

To prove the positivity of ips for e small, we first establish that limsupg_^g l-B^I 
is finite, from which we deduce that the constants A^ are positive; in fact, 
we will find that A^ > 1 — 45. Let us suppose already that 1 — 55 > and 
Ai(£) < 2A* - 5 for all e e £2. Since ip^ e Cx,{R), for alH < 0, £ e £^2 \ {0}, 
we get 

\A,cxp{Xi{e)t) + B,cxp{X^{e)t)\ < \ij,{t)\ + \xs{t)\ < C^cxpiXj), 
where C5 — Ci + C4. In particular, taking t — and t — —1, we obtain 
\A,+B,\ < C5, \A,+B,exp(Xi{s)-X^(e))\ < C5exp(Ai(£)-A*) < C5exp(A), 
hence \B^\{1 - exp(Ai(£) - X^{s))) < C^il + exp(A)) := Cq, for £ G £2 \ {0}. 

Noting that exp(Ai(£) — Aoo(^)) ^ as £ ^ 0, we deduce that there is 
£s — (—£3, £3) C S2 such that l-Bel < ^Cq for £ e £3 \ {0}, so that 

\B,exp{^oo{e)t)\ < 2C6exp(Aoo(£)t) < 2Cg &q){5t) &q){{2X, - 5)t) , 

for t < 0, e e SsX {0}. Set y^{t) = B^exp{Xoo{e)t) + Xs{t). By Lemma 8, 
we have ■^(t) = exp(At) + z{t) with z{t) — 0(exp(2A*t)) at t — —00. Since 
limt_>_oo Cq exp{5t) — 0, we now conclude that there is sq — so{5) < such 
that for t < So and < |£| < £3 we have 

\ye{t)\ < 5exp((2A, - 6)t), \ye{t) - z{t)\ < (5exp((2A* - 6)t). 

On the other hand, for So = Sexp{{X — A*)so) = 6exp{6so), there exists 
£4 = £4(5) e (0,£3] such that, for |£| < £4, we have \ipsit) —ip{t)\ < 5oexp{X^t). 
Taking t — So we obtain for < |£| < £4 

|^£exp(Ai(£)so)-exp(Aso)| < \ipe{so)-'ip{so)\ + \ye{so)-z{so)\ < 25 exp(Aso), 

hence ^e(so) = exp(Ai(£)so) + 2/e(so) > exp(Aso) - (l^lg exp(Ai(£)so) - 
exp(Aso)| + \ye{so)\) > (1 ~ 35)exp(Aso). Therefore, for all < |£| < £4, 

Ae > -Veiso) exp(-Ai(£)so) + (1 - 35) exp((A - Ai(£))so) > (1 - 45) > 0. 
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Thus, for < |£| < £4 and t < sq we get ipeit) > exp{Xi{e)t)[{l - A5) - 5)] > 0. 
Since linij^o "^e = uniformly on M and ip is bounded from below by a positive 
constant on [sq, 00), we conclude that ip^ is positive on M, for all s small. 

Finally, for every fixed £ G £2 \ {0}, we have that g{x) = px + q{x)x'^, ipeit) — 
exp(Ai(£)t) + exp(2A^,t), where q G C[0, +cxo) and is bounded on 
(—00,0]. Hence, '0£(t) = ^£ exp(Ai(£)i) + 0(exp(2A*i)) at —00 and 

g{i\)^{i - h)) = A,pexp(-Xi(e)h) exp(Xi{e)t) + Cs(t) exp(2A,i), e e {0}, 

where Ce{t) is bounded: \ce{t)\ < co{s), t <0. Differentiating (20), we obtain 



\ —00 

H — J e 5? g{i/je{s - h))ds 

t 

\ —00 

+ ^J;^ j e ^ exp(Ai(£)s)ds +0(exp(2A,t)) = 

= As\i{e) exp(Ai(£)i) + 0(exp(2A*i)), t -00. 

Hence, ipe,-ip'e ^ C'ai(£)(M) so that ^'^ G C2Ai(£)(M) in (25). Therefore, in 
view of [24, Proposition 7.1] and the inequahty A < Xi{e), e G 82 \ {0}, 
we get from (25) that i^^it) = exp(Ai(£)i) + 0(exp(2Ai)) and = 
A^Xi{e) exp(Ai(£)i) + 0(exp(2Ai)), t -00. □ 




4 Non-monotonicity of travelling waves 



As it was noticed in [1,9,10,12,20,29], various investigators have studied nu- 
merically the case of a large delay in the Nicholson's blowflies equation, and 
noted a loss of monotonicity of the wave front as the delay increases, "with 
the front developing a prominent hump " whose height "is bounded above by a 
bound that does not depend on the delay", see [1, p. 308] from which the above 
citation was taken. It is not difficult to explain the second phenomenon, since 
at every point of local maximum a of ijj{t, c) wc have ip'{(y, c) = 0, ip"{(J, c) < 
so that ipi^a, c) < g{ip{o' — h, c)) < max^>o5'(x). Here wc explain also the first 
phenomenon, easily getting the oscillation of the travelling waves about K as 
t — > +00 stated in Theorem 1 from the next two lemmas. 
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Lemma 15 Let g'{K) <0,h> and \g'{K)\he'^'^^ > 1. Then the equation 



has no negative real roots, for all sufficiently small e. Moreover, if the equilib- 
rium K of (4) is hyperbolic, then, for all small e, there are no roots of (28) 
on the imaginary axis. 

PROOF. Set /^e{z) = £^z^ -z-l + g'{K) exp{-zh). We first prove that the 
lemma is vahd for £ = (see also [14]). Let Zq be the maximum point of Ao{z) 
on M, i.e., £ is such that Aq(2;o) = —l — hg'{K) exp(— ^o/i) = 0. Note that 
Ao(— c)o) = —oo and Aq{z) < for 2; > 0. If there is a negative zero of Ao{z), 
then -20 < and Ao(^o) = —zq — 1 — l/h > 0, implying that = Aq(^o) > 
— 1 + \g'{K)\he^'^^, which contradicts the hypothesis \g'{K)\he'^'^^ > 1. 

Uh\g'{K)\ > 1, then A^(^) > 2ze^-l + e-''' > z{2e^-h) > for all ^ < and 
< h/2, hence As{z) < for ^ < 0. Now, let h\g'{K)\ < 1, so that zq < 0. 
For > small, by the implicit function theorem we conclude that there is a 
negative root z{e) of the equation A^(z) = with z{0) — Zq] moreover, z{e) is 
the absolute maximum point of ^ '^^{z) on (— cxo, 0]. Since S{e) := As{z{e)) 
depends continuously on e and 5(0) < 0, for £ > small we have ^^{z) < 
for all z <0. 

We now prove that (28) has no roots on the imaginary axis. First, notice that 
\Ae{ib)\ > \g'{K)\ > if 6 > 2\g'{K)\. For e = 0, Eq. (28) does not have 
roots on the imaginary axis, therefore |Ao(?6)| > 0, |6| < 2\g'{K)\. Hence, 
|A£(i6)| > for all e small and |6| < 2\g'{K)\, which implies the hyperbolicity 
of Eq. (28). □ 

The next lemma can be considered as an extension of the linearized oscillation 
theorem from [15] to the second order delay differential equation 



Lemma 16 Assume (H) and that g'{K)he^'^^ < —1. For small e > 0, set 

{Dx){t) = e'^x"{t) —x'{t) —x{t) + g'[K)x(t — h). Then every non constant and 
bounded solution x : M ^ M 0/ (29) such that a;(+oo) = K oscillates about K. 

PROOF. Consider some non-constant solution x : R ^ R of (29) such that 
a;(+oo) = K. If we suppose for a moment that, for some 77 G M, it holds 
x{s) = K identically for all s > rj, then we obtain easily that x{s) = K for all 
s G [?7 — /i, r;]. Hence x should be a constant solution, in contradiction with our 
initial assumption. Therefore either a{t) = x{t) — K oscillates about zero or 



e^z^ -Z-1 + g\K) exp(-^/i) = 



(28) 



e^x"{t) - x'{t) - x{t) + g{x{t -h))^0, te M. 



(29) 
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is eventually non-constant and non-negative, or non-positive. In order to get 
a contradiction, assume that a is not oscillatory. Notice that a satisfies the 
following linear asymptotically autonomous delay difi^erential equation 

- (7'{t) - a(t) + -f{t)a(t - /i) = 0, teR, 7(+oo) = g'(K) < 0, (30) 

where 7(t) = g'{K) + co{t) and co{t) = g\K + e{t)a{t - h)) - g'{K) for 
some 9{t) e (0, 1) given by the mean value theorem. Since x{t) is bounded 
on M and x{t) ^ K as t +oo, we can use the integral representations 
(20) and (27) to prove that limt^^^ x' (t) = 0. From Lemma 15, it follows 
that the equation {Dx){t) = is hyperbolic, hence the equilibrium {K,0) of 
the system x'{t) — v{t),s'^v'{t) — v{t) — x{t) + g{x{t — h)) = is hyperbolic 
for all sufficiently small e. Thus the trajectory of x{t) belongs to the stable 
manifold of the hyperbolic equilibrium K of (29), so that we can find a > 
such that a{t) = 0(e~"*), cr'(t) = 0(e~"*) at t = +cxo. Therefore we have 
co(t) = 0{a{t - h)) = 0{e-"-^) at +oo. From [24, Proposition 7.2] (see also 
[18, Proposition 2.2]), we conclude that: (i) either there are h> a, (5 > and 
u{t) a nontrivial eigensolution of the limiting equation 

e^u"{t) - u'{t) - u{t) + g\K)u{t - h) = (31) 

on the generalized cigcnspace associated with the (nonempty) set A of eigen- 
values with 3?e A = —6, such that ait) = u(t) + 0(exp(— (6 + 6)t), t — +cxo; 
{ii) or a{t) decays superexponentially at +oo. However, this latter condition 
is not possible: as it was estabhshed in [19, Lemma 3.1.1], if '^{+00) 7^ then 
every eventually nontrivial and nonnegative solution of (30) does not decay 
superexponentially (see also [18, Lemma A.l] for the case 7(-|-cxo) > 0). On the 
other hand, from Lemma 15 we know that there are no real negative eigenval- 
ues of (31): hence Q'm A 7^ for all A G A. From [18, Lemma 2.3], we conclude 
that a{t) is oscillatory. □ 

Finally, we observe that due to the exponential stability of the positive steady 
state, which implies fast convergence, numerical heteroclinic solutions '^(i, c) 
exhibit only one or two well pronounced humps, see [9, Fig. 2]. 
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